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ABSTRACT 

Recent X-ray observations have been used to demonstrate that the cluster gas mass - temperature 
relation is steeper than theoretical self-similar predictions drawn from numerical simulations that consider 
the evolution of the cluster gas through the effects of gravity and shock heating alone. One possible 
explanation for this is that the gas mass fraction is not constant across clusters of different temperature, 
as is usually assumed. Observationally, however, there is no compelling evidence for gas mass fraction 
variation, especially in the case of hot clusters. Seeking an alternative physical explanation for the 
observed trends, we investigate the role in the cluster gas mass - temperature relation of the preheating of 
the intracluster medium by some arbitrary source for clusters with emission- weighted mean temperatures 
of Tx ^ 3 keV. Making use of the physically-motivated, analytic models developed in 2002 by Babul and 
coworkers, we find that preheating does, indeed, lead to a steeper relation. This is in agreement with 
previous theoretical studies on the relation. However, in apparent conflict with these studies, we argue 
that a "high" level of entropy injection is required to match observations. In particular, an entropy floor 
of > 300 keV cm^ is required. We also present a new test, namely, the study of the relation within 
different fixed radii. This allows one to indirectly probe the density profiles of clusters, since it samples 
different fractions of the virial radius for clusters of different temperature. This test also confirms that 
a high level of preheating is required to match observations. 

Subject headings: cosmology: theory — galaxies: clusters: general — X-rays: galaxies: clusters 



1. INTRODUCTION 

Analytic models and numerical simulations of clusters 
of galaxies have been used to predict the existence of scal- 
ing relations between various observable quantities, such 
as the well-known luminosity (ix) - temperature {Tx) 
and mass (M) - temperature relations, where Lx oc T^ 
and M oc T_^'^, respectively. However, it is now fairly well 
established that X-ray properties of clusters do not scale 
in such a fashion. Most notable of these is the Lx — Tx 
relationship, which is observed to be much steeper than 
predicted, Lx oc r|-6-3.o (^g_g^^ Markevitch 1998; Allen & 
Fabian 1998; Arnaud & Evrard 1999). 

Considerable effort has recently been directed towards 
explaining why the observed relations deviate from their 
predicted scalings (e.g., Tozzi & Norman 2001; Dave et al. 
2001; Babul et al. 2002, hereafter BBLP02). In particular, 
it is the Lx — Tx relation that has grabbed most of the 
spotlight because there is a wealth of published observa- 
tional studies on the luminosities and temperatures of clus- 
ters with which to compare models and simulations. How- 
ever, another important scaling relation is the cluster gas 
mass (Mgas) - Tx relation. Neumann & Arnaud (2001) 
have suggested that a deviation from the self-similar scal- 
ing of Mgas oc T^^ might "explain" the observed deviation 
in the Lx —Tx relation. Indeed, a number of observational 
studies have indicated that the relation is much steeper, 
with Mgas oc (vikhhnin et al. 1999; Mohr et 

al. 1999, hereafter MME99; Neumann & Arnaud 2001). 
If the gas density profile is roughly self-similar, this does 
lead to consistency with the observed Lx — Tx relation. 
However, we still need a physical explanation for why the 
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relationship between a cluster's gas mass and its temper- 
ature deviates from its self-similar scaling. 

Expressing the total gas mass within the cluster as 
Mgas = fgasM, a Steepening of the Mgas - Tx relation 
can be interpreted as a dependence of fgas on cluster 
mass. That is, if M oc T^^ , as suggested by the self-similar 
model, then the observed Mgas — Tx relation implies that 
fgas oc T^'^""'^. A varying gas mass fraction is expected 
if the efhciency of galaxy formation varies systematically 
across clusters of different mass. Observational support for 
this has been claimed recently by Bryan (2000). However, 
this is still controversial, and there is no compelling evi- 
dence for a variation of fgas with cluster temperature (but 
see Arnaud & Evrard 1999; MME99). This is especially 
true for the systems that we are specifically interested in: 
hot clusters with Tx ^ 3 keV. This is apparent, for exam- 
ple, in Figure 1 (top) of Balogh et al. (2001), who carry 
out an accounting of stars and gas to estimate the fraction 
of coohng baryons in clusters. Moreover, Roussel, Sadat, 
& Blanchard (2000) have carried out a careful analysis of 
group and cluster X-ray data to estimate fgas directly and 
have found no trends. More recently, Grego et al. (2001) 
have analysed Sunyaev-Zeldovich effect observations of 18 
hot clusters and have also found no correlations between 
a hot cluster's gas mass fraction and its temperature. Fi- 
nally, observational studies of the total cluster mass (M) 
- temperature relation have indicated that M oc y^-S~2.o 
(Horner et al. 1999; Ettori & Fabian 1999; Nevalainen et 
al. 2000; Finoguenov et al. 2001), which, given the ob- 
served Mgas ~ Tx relation, is consistent with fgas being 
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constant. 

Theoretically, it is only now becoming possible to reli- 
ably investigate the dependence of fgas on temperature 
with the inclusion of radiative cooling, star formation, 
feedback, and other relevant processes in numerical simu- 
lations (e.g., Lewis et al. 2000; Pearce et al. 2000; Muan- 
wong et al. 2001; Dave et al. 2001). As of yet, how- 
ever, there is little agreement in the approaches adopted 
to model these processes and prevent the so-called cool- 
ing crisis (compare, for example, the findings of Lewis ct 
al. 2000 with those of Pearce et al. 2000). This is not 
surprising. As discussed in detail by Balogh et al. (2001), 
attempting to model the effects of cooling across the wide 
range of halo masses found in clusters is inherently very 
difficult. The addition of "sub-grid" processes, such as 
star formation and feedback, further complicates matters. 
Thus, the effects that these additional physical processes 
have on the gas mass fraction of clusters will not be fully 
realized until such issues are resolved. 

In this paper, however, we show that the observed vari- 
ation of the Mgas — Tx relation(s) arises quite naturally 
within the class of models that invoke preheating of the 
intraclustcr medium during the early stages of cluster for- 
mation. In these models, fgas is constant on cluster scales 
{'^x ^ 3 keV), and the self-similarity is instead broken 
by an entropy floor generated by early non-gravitational 
heating events. Preheating has previously been shown to 
bring consistency between a number of other observed 
and predicted scaling relations for groups and clusters 
(e.g., BBLP02), and therefore one might expect that the 
Mgas — Tx relation should also be modified. 

The preheating model was originally put forward by 
Kaiser (1991) and has subsequently been investigated by 
a number of authors (e.g., Evrard & Henry 1991, Bower 
1997, Cavaliere et al. 1997; 1998; 1999; Balogh et al. 1999, 
Wu, Fabian, & Nulsen 2000; Locwenstein 2000, Tozzi & 
Norman 2001; Borgani et al. 2001; Thomas et al. 2002; 
BBLP02). If the ICM is injected with enough thermal en- 
ergy, the hot X-ray emitting gas will become decoupled 
from the dark halo potential and break the self-similar 
scaling relations. The best estimates suggest that a sub- 
stantial amount of energy (~ 1 keV per particle) is required 
to reproduce the observed relations (mainly the Lx — Tx 
relation). It is not yet known what source(s) could in- 
ject such a large amount of energy into the ICM. Both 
galactic winds (driven by supernovae) and ejecta from ac- 
tive galactic nuclei have been proposed, but because of 
the complexity of the physics, the exact details have yet 
to be worked out. For an in-depth discussion of potential 
sources of preheating and of alternative possibilities for 
reproducing the observed relations we refer the reader to 
BBLP02. 

In this paper, we adopt the physically motivated, ana- 
lytic model developed in BBLP02 to explore the impact 
of cluster preheating on the Mgas — Tx relation. In com- 
parison with the Lx — Tx and M — Tx relations, it has 
drawn very little attention by theoretical studies. The only 
studies to have examined the effects of entropy injection on 
the Mgas — Tx relation to date are Locwenstein (2000) and 
Bialek et al. (2001). To be specific, Locwenstein (2000) 
considered models where the entropy injection occurs at 
the centers of groups and clusters, after the latter have 



formed whereas Bialek et al (2001), like BBLP02, investi- 
gated preheated models. For reasons that will be described 
below (in §5), we believe our work greatly improves upon 
both of these studies. The prevailing apathy by theorists is 
perhaps due in part to a near absence of published obser- 
vational studies on the gas masses of clusters. However, in 
light of the recent important observations discussed above 
(e.g., Vikhlinin et al. 1999; MME99; Neumann & Ar- 
naud 2001) and the new influx of high resolution data from 
the Chandra and XMM-Newton X-ray satellites, which will 
likely provide even tighter constraints, we believe that a 
thorough examination of the Mgas — Tx relation is timely. 

The models we consider below were developed in a flat 
A-CDM cosmology with ilm = 0.3, h = 0.75, and a nucle- 
osynthesis value hb = 0.019/1-2 (Buries & Tytler 1998). 
They are computed for a number of different preheating 
levels, corresponding to entropy constants of Kq = 100, 
200, 300, & 427 keV cm^. These span the range required 
to match the observed Lx — Tx relations of groups and 
hot clusters (e.g., Ponman et al. 1999; Lloyd-Davies et al. 
2000; Tozzi & Norman 2001; BBLP02). For the purposes 
of comparison, we also implement an "isothermal" model 
(see Section 2.3 in BBLP02), which mimics the self-similar 
result deduced from numerical simulations (e.g., Evrard et 
al. 1996). 

2. CLUSTER MODELS 

Since an in-depth discussion of the preheated cluster 
models can be found in BBLP02, we present only a brief 
description of the models here. 

The preheated models can be summarized as follows: 
the dominant dark matter component, which is unaffected 
by the energy injection, collapses and virializes to form 
bound halos. The distribution of the dark matter in such 
halos is assumed to be the same as found in recent ultra- 
high resolution numerical simulations (Moore et al. 1998; 
Klypin et al. 1999; Lewis et al. 2000) and is described by 

Pdm{r) = Pdmfi ^1 + (1) 

where n = 1.4, pdm.o is the profile normalization, and Ts is 
the scale radius. While the dark component is unaffected 
by energy injection, the collapse of the baryonic compo- 
nent is hindered by the pressure forces induced by preheat- 
ing. If the maximum infall velocity due purely to gravity 
of the dark halo is subsonic, the flow will be strongly af- 
fected by the pressure and it will not undergo accretion 
shocks. It is assumed that the baryons will accumulate 
onto the halos isentropically at the adiabatic Bondi ac- 
cretion rate (as described in Balogh et al. 1999). This 
treatment, however, is only appropriate for low mass ha- 
los. If the gravity of the dark halos is strong enough (as 
it is expected to be in the hot clusters being considered 
here) so that the maximum infall velocity is transonic or 
supersonic, the gas will experience an additional (gener- 
ally dominant) entropy increase due to accretion shocks. 
In order to trace the shock history of the g detailed 
knowledge of the merger history of the cluster /group is 
required but is not considered by BBLP02. Instead, it 
is assumed that at some earlier time the most massive 
cluster progenitor will have had a mass low enough such 
that shocks were negligible in its formation, similar to the 
low mass halos discussed above. This progenitor forms an 
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iscntropic core of radius Tc at the cluster center. The en- 
tropy of gas outside of the core, however, will be affected 
by shocks. Recent high resolution numerical simulations 
suggest that the entropy profile for gas outside this core 
can be adequately represented by a simple analytic ex- 
pression given by In K(r) = IhKq + aln(r/rc) (Lewis et 
al. 2000), where a ~ 1.1 for the massive, hot clusters 
{Tx > 3 keV) of interest here (Tozzi & Norman 2001; 
BBLP02). It should be noted that in the case of these 
massive systems, the accretion of gas is limited by grav- 
itational infall, and hence they accrete their full compli- 
ment of baryons [i.e., Mgas = {^b/^m)M]. It is assumed 
that the mass of baryons locked up in stars is negligible 
(as suggested by, for example, Rousscl et al. 2000; Balogh 
et al. 2001). 

Following this prescription and specifying the parame- 
ters Tc, Pgas{rc), and a (as discussed in BBLP02) com- 
pletely determines the models. Under all conditions, the 
gas is assumed to be in hydrostatic equilibrium within the 
dark halo potential. The effects of radiative cooling are 
neglected by these models. 



3. HOW PREHEATING AFFECTS THE mgas - tx RELATION 

Preheating will affect the Mgas — Tx relation in two 
ways: (1) by altering the temperature profile and increas- 
ing the emission- weighted gas temperature of a cluster and 
(2) by altering the gas density profile and reducing the gas 
mass in the cluster core. We are interested in the strength 
of these effects and whether or not they can be distin- 
guished by current or future observational data. First, we 
consider the effect of preheating on the temperature of a 
cluster. 

Figure 1 is a plot of Tx as a function of entropy floor 
level (i.e., Kq) for three clusters of different total masses. 
The thin line represents a cluster with M « 5.6 x 1O^*M0, 

the next thickest line represents a cluster with M « 
IO^^Mq, and the thickest line represents a cluster with 
M w 1.8 X IO^^Mq. As expected, the gas temperature of 
a cluster increases as the level of preheating is increased. 
On average, an increase in Tx of about 1 keV (10 - 25%) 
occurs when a cluster is preheated to the level of Kq > 300 
keV cm2 (over the range 3 keV <Tx <10 keV). This ef- 
fect will primarily manifest itself as a normalization shift 
in the Mgas — Tx relation. 

Figure 2 presents the dimensionlcss gas density profile 
of a cluster with Tx = 4 keV (left panel) and a cluster 
with Tx = 8 keV (right panel) as a function of the level 
of preheating. The dot-dashed line is the self-similar re- 
sult (i.e., isothermal model of BBLP02). The long-dashed, 
short-dashed, dotted, and solid lines represent the pre- 
heated models of BBLP02 with Kq = 100, 200, 300, and 
427 keV cm^, respectively. Preheating reduces the gas 
density, and therefore the gas mass, in the central regions 
of a cluster. In §4, we investigate the Mgas — Tx rela- 
tion within three different radii: r = 0.25h~^ Mpc and 
0.50/i~^ Mpc and r^oo (the radius within which the mean 
dark matter mass density is 500 times the mean critical 
density Pcrit at z = 0). These radii arc indicated in Figure 
2 by the open squares, pentagons, and triangles, respec- 
tively. Clearly, the effect on the Mgas — Tx relation will be 
strongest when Mgas is evaluated within r = 0.25h~^ Mpc. 
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Fig.l. Effect of preheating on a cluster's temperature. Tfie thin 
line, the next thickest line, and the thickest line represent clusters 
with M f» 5.6 X 10", lO^^, and 1.8 x W^^Mq, respectively. The 
solid squares indicate the discrete points where the model was actu- 
ally evaluated. 

Furthermore, because r ~ 0.25h~^ Mpc is a fixed ra- 
dius that samples different fractions of the virial ra- 
dius {Rfiaio) for clusters of different temperature, the 
(fractional) reduction in gas mass within that radius 
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Fig. 2. Effect of preheating on a cluster's gas density profile. Left: 
Cluster with Tx = 4 keV. Right: Cluster with Tx = S keV. The dot- 
dashed line is the self-similaj: result. The long-dashed, short-dashed, 
dotted, and solid lines represent the preheated models of BBLP with 
Kq = 100, 200, 300, and 427 keV cm^, respectively. The squares, 
pentagons, and triangles indicate the radii r = 0.25h~^ Mpc and 
0.50h~^ Mpc and rsoo, respectively, for each of the models. 
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will be largest for the lowest temperature systems. This 
will lead to both a normalization shift and a steepening 
of the Mgas — Tx relation. To illustrate how strong the 
effect is, we examine the reduction of gas mass within 
r — Q.2bh~^ Mpc, Q.hOh~^ Mpc and rsoo for a low mass 
cluster with a total mass of M f» 5.6 x lO^^M© and for a 
high mass cluster with M Ri 1.8 x lO^^M©. We find that 
when the low mass cluster has undergone preheating at the 
level of fsTo = 100 keV cm^ it has « 32% more mass in gas 
within r = 0.25/i^"'^ Mpc than when the same cluster has 
undergone preheating at the level of Kq = 427 keV cm^. 
Using the same test on the M w 1.8 x lO^^M© cluster, 
however, yields a difference of only 22%. When we probe 
the larger radius r = 0.50/i~^ Mpc, we find the effect is 
less pronounced (as expected). The difference in Mgas be- 
tween the Kq = 100 keV cm^ model and Kq = 427 keV 
cm^ model is 8% for the low mass cluster and 7% for the 
high mass cluster. Finally, when the gas mass is evaluated 
within rsoo, the difference is 4% for the low mass cluster 
as opposed to 2% for the high mass cluster. 

In summary, preheating will significantly affect the 
Mgas — Tx relation by increasing the emission-weighted 
gas temperature of clusters. Whether or not the relation 
is also affected by the reduction of gas mass in the cores 
of clusters depends on within which radius Mgas is evalu- 
ated and what temperature regime is being probed. The 
effect will be strongest for low temperature systems and 
when Mgas is probed within small radii (e.g., r = Q.25h~^ 
Mpc). An evaluation of the Mgas — Tx relation within 
large radii, such as rsoo, however, probes the integrated 
properties of a cluster and will be sensitive only to the 
temperature shift. 

In the next section, we compare the results of the 
BBLP02 preheated models to genuine observational data. 
As we show below, only models with Kq > 300 keV cm^ 
are consistent with the data. 



4. RESULTS 

In Figure 3 wc present the Mgas — Tx relation as pre- 
dicted by the BBLP02 preheated models within rsoo. The 
radius rsoo is typically comparable in size to the observed 
radius of a cluster and represents the boundary between 
the inner, virialized region and recently accreted, still set- 
tling outer region of a cluster (Evrard et al. 1996). Thus, 
as already mentioned, the Mgas{r 500) — Tx relation can 
be regarded as a probe of the integrated properties of a 
cluster and can be directly compared with the self-similar 
result oiMgas ocT^-^. 

In Figures 4 and 5 wc present the Mgas — Tx rela- 
tion as predicted by the BBLP02 preheated models within 
the fixed radii r = 0.25/1"^ Mpc and r = 0.50/1"^ Mpc, 
respectively. As mentioned above, the determination of 
the Mgas — Tx relation within some fixed radius, such as 
r = 0.25h~^ Mpc or r = 0.50h~^ Mpc, can be used as 
an indirect probe of the gas density profiles of clusters be- 
cause it samples different fractions of the virial radius for 
clusters of different temperature. For the purposes of clar- 
ity, we discuss the Mgas — Tx relation at these three radii 
separately. 
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Fig. 3. Comparison of Mgasi^soo) ~ Tx relations. Squares rep- 
resent the gas mass determinations of Mohr ct al. (1999) within 
rsoo. The dot-short-dashed, long-dashed, short-dashed, dotted, and 
solid lines represent the self-similar model and preheated models of 
BBLP02 with Kq = 100, 200, 300, and 427 kcV cm^, respectively. 
The dot-long-dashed (thick) line represents the best-fit heated model 
of Loewenstein (2000). 
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Fig. 4. Comparison of Mgas (r = 0.25fe~i Mpc)-Tx relations. The 
solid triangles and pentagons represent the gas mass determinations 
Peres et al. (1998) and White et al. (1997) within r = 0.25/1-^ 
Mpc, respectively. The dot-dashed line is the self-similar result. The 
long-dashed, short-dashed, dotted, and solid lines represent the pre- 
heated models of BBLP02 with Kq = 100, 200, 300, and 427 keV 
cm^, respectively. 
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Fig. 5. Comparison of Mgasir = 0.50h^^ Mpc)-Tx relations. The 
solid pentagons represent the gas mass determinations of White et 
al. (1997) within r = O-bOh'^ Mpc. The dot-dashed line is the 
self-similar result. The long-dashed, short-dashed, dotted, and solid 
lines represent the preheated models of BBLP02 with Kq = 100, 200, 
300, and 427 keV cm-^, respectively. 

4.1. Test 1: Mgasirzoo) - Tx 

The solid squares in Figure 3 represent the gas mass de- 
terminations of MME99 within r.soo using surface bright- 
ness profile fitting (with isothermal /3 models) of ROSAT 
Position Sensitive Proportional Counter data and mean 
emission-weighted temperatures from the literature, for 
clusters with Tx > 3 ke'V and whose error bars are 1 keV 
or smaller. We compare this with the self-similar result 
represented by the "isothermal" model of BBLP02 (dot- 
short-dashed line). Finally, the long-dashed, short-dashed, 
dotted, and solid lines represent the preheated models of 
BBLP02 with Kq = 100, 200, 300, and 427 keV cm^, re- 
spectively. The thick dot-long-dashcd line represents the 
predictions of the best-fit heated model of Loewenstein 
(2000). This model is discussed further in §5.1. 

It is readily apparent that only the preheated models of 
BBLP02 with Kq > 200 keV cm^ have a reasonable chance 
of being consistent with the data of MME99. The normal- 
ization clearly indicates that the observed gas temperature 
of clusters with a given gas mass is hotter than predicted 
by models with entropy floors of Kq < 100 keV cm^. We 
note that this discrepancy may be remedied by assuming 
a smaller value of Qb/Qm- However, a similar offset, in 
the same sense, is seen in the correlation with total dark 
matter mass and gas temperature (Horner et al. 1999; 
Nevalainen et al. 2000; Finoguenov et al. 2001). This will 
not be reconciled by lowering Qb/^m- The reason why the 
preheated models with Kq > 200 keV cm^ are better able 
to match the normalization of the observational data than 
models with Kq < 100 keV cm^ is, as mentioned above, 
because an increase in the amount of preheating directly 

^ Our best fit differs slightly from MME99's best fit to their own data 
used only clusters with Tx > 3 keV and whose error bars are 1 keV or 



leads to an increase in the emission-weighted gas temper- 
ature. 

We have attempted to quantify how well (or poorly) 
the preheated and self-similar models match the obser- 
vational data. We have fit both the theoretical results 
and the observational data with simple linear models of 
the form logMg^g = mlogTx + b over the range 3 keV 

< Tx < 10 keV. For the theoretical results, we calcu- 
late the best-fit slope and intercept using the ordinary 
least squares (OLS) test. We stress that the results of 
these fits, which are presented in Table 1, are only valid 
for clusters with Tx > 3 keV. At lower temperatures, the 
role of preheating becomes much more important [as Mgas 
becomes less than {fti,/flm)M] and as a result, the rela- 
tions steepen dramatically. For example, the preheated 
model with Kq = 427 keV cm^ is well approximated by 
a power-law with Mgas oc Tj^^^ over the range 3 keV 

< Tx ^10 keV but is signiflcantly steeper over the range 
1 keV < Tx < 3 keV with Mgas « T^-^^. Thus, it is 
absolutely essential that comparisons between theoretical 
models and observations are done over the same range in 
temperatures. 

To fit the observational data of MME99, we have used 
a linear regression technique that takes into account mea- 
surement errors in both coordinates as well as intrinsic 
scatter (the BCES test of Akritas & Bershady 1996). As 
a consistency check, we have also employed 10,000 Monte 
Carlo bootstrap simulations. No significant deviations be- 
tween the two tests were found. The results of the linear 
regression fits to the observational data are also presented 
in Table 1. 

For all 38 clusters taken from MME99, we derive a best 
fit that is inconsistent with the results of all of the theo- 
retical models considered at greater than the 90% confi- 
dence level^. However, as is apparent from Figure 3, the 
slope and intercept of the best-fit line are sure to be heav- 
ily dependent upon the two low temperature clusters with 
the lowest measured gas masses (and gas mass fractions): 
the Hya I cluster (Abell 1060) and the Cen cluster (Abell 
3526). A number of other studies (both optical and X- 
ray) have also identified very unusual properties in both 
clusters. For example, Fitchett & Merritt (1988) were un- 
able to fit a spherical equilibrium model to the kinematics 
of galaxies in the core of Hya I. They suggest that sub- 
structure is present and is likely why Hya I does not lie 
along the Lx — o" relation for galaxy clusters. More re- 
cently, Furusho et al. (2001) have found that the metal 
abundance distribution implies that the gas in Hya I is 
well-mixed (i.e., it does not contain an obvious metallicity 
gradient), suggesting that a major merger event may have 
occurred sometime after the enrichment of the ICM. Mea- 
surements of the bulk motions of the intracluster gas in 
the Cen cluster (through Doppler shifting of X-ray spec- 
tral lines) reveal strange gas velocity gradients indicative 
of a large merger event in the not too distant past (Dupke 
& Bregman 2001). This picture has also been supported 
by Furusho et al. (2001), who found large temperature 
variations across the cluster's surface. Thus, neither Hya 
I nor Cen can be regarded as typical "relaxed" clusters 
and are probably not representative of the majority of low 

because we implemented a different selection criteria. Namely, we have 
smaller. 
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temperature systems. 

One way to ameliorate the impact of the two clusters 
would be to increase the number of systems of this tem- 
perature. However, there are very few published gas mass 
estimates of cool clusters and groups within r^ao- X-ray 
emission from groups is usually only detected out to a 
small fraction of this radius. The one study that does 
present group gas masses for a radius at fixed overdensity, 
Roussel et al. (2000), does so for r2oo and is not directly 
comparable to the results presented in Figure 3. Also, 
in that study, gas masses were determined by extrapolat- 
ing the surface brightnesses far outside the limiting radius 
for which X-ray emission was actually detected. This can 
lead to biases in determining group/cluster properties (see 
Mulchaey 2000: Balogh et al. 2001). 

In recognition of the above, we have tried removing Hya 
I and Cen from the sample and fitting the remaining 36 
clusters using the same procedure. We find that the pre- 
heated model with Kq = 427 keV cm^ is then consistent 
with the data at the 90% level. The Kq = 200 and 300 kcV 
cm^ models are marginally inconsistent with the MME99 
data. The isothermal model is ruled out at > 99% con- 
fidence irrespective of whether these clusters are dropped 
or not. 

The other two observational studies that have investi- 
gated the Mgas — Tx relation, Neumann & Arnaud (2001) 
and Vikhlinin et al. (1999), unfortunately did not present 
gas mass determinations within rsoo for individual clus- 
ters in a table or graphically. They did, however, present 
their best-fit values for the slope of the relation. These 
were deduced from samples of clusters that have tempera- 
tures spanning roughly the same range as that considered 
in Figure 3. The best-fit slopes of the preheated mod- 
els are shallower than the best-fit claimed by Neumann & 
Arnaud (2001) of Mgas oc T^-^^ for a sample of 15 hot clus- 
ters. However, an estimate of the uncertainty on this result 
was not reported; thus we are unable to say whether this 
result is inconsistent with the predictions of the preheated 
models. The predicted slopes of all four preheated models 
studied here are in excellent agreement with the findings of 
Vikhlinin et al. (1999), who report Mgas oc T^-7i±o i3 f^^. 
their sample of 39 clusters. We also note that the results 
of Neumann k Arnaud (2001) and Vikhlinin et al. (1999) 
differ significantly from the predictions of the self-similar 
model. 

In summary, we find that the class of models that invoke 
preheating are much better able to match the observed 
Mgas {^500) — Tx of hot clusters than that of the isother- 
mal self-similar model, which is ruled out with a high level 
of confidence. A careful analysis of the MME99 data also 
suggests that only those models that invoke a "high" level 
of energy injection (i.e., Kq > 300 keV cm^) are able to 
match observations. 

4.2. Test 2: Mgas(r = 0.25/1"! Mpc) - Tx 

The solid triangles and pentagons in Figure 4 represent 
the gas mass determinations within r = 0.25/i~-^ Mpc of 
Peres et al. (1998) and White et al. (1997), respectively. 
These data were obtained using surface brightness profile 
fitting of ROSAT data (Peres et al. 1998) and Einstein 
data (White et al. 1997) and emission-weighted tempera- 
tures from the literature, for clusters with Tx > 3 keV and 



whose error bars are 1 keV or smaller. The predictions 
of the isothermal self-similar model are represented by 
the dot-dashed line. Once again, the long-dashed, short- 
dashed, dotted, and solid lines represent the preheated 
models of BBLP02 with Kq = 100, 200, 300, and 427 keV 
cm^, respectively. 

In spite of the scatter, it is apparent that only those pre- 
heated models with entropy floors of Kq > 300 keV cm^ 
are consistent with the 57 clusters plotted in Figure 4. As 
with the Mgasir^oa) —Tx relation, the normalization of the 
self-similar model and the preheated model with Kq = 100 
keV cm^ suggests that ICM is observed to be much hotter 
than predicted by either of these models. Fitting both the 
theoretical predictions and observational data in a manner 
identical to that presented in the previous subsection, we 
find that only the preheated models with Kq > 300 keV 
cm^ have both slopes and intercepts that are consistent 
with the observational data (see Table 2) . On the basis of 
normalization (intercept), the self-similar model is ruled 
out with greater than 99% confidence. 

In §3 we briefly discussed the potential of the gas density 
profile to affect the Mgas{r = 0.25h~^ Mpc) - Tx relation. 
This effect is obvious in Figure 4, with mild breaks at 
Tx « 10 keV for the Kq = 427 keV cm^ model and at 
Tx « 5 keV for the Kq = 300 keV cm^ model. However, 
with the large scatter obscuring any potential breaks in 
the Mgas — Tx relation, all we can conclude is that the 
data are consistent with predicted profiles of the BBLP02 
preheated models with Kq > 300 keV cm^. 

The exact nature of the scatter in Figure 4 is unclear. 
While some of the scatter is likely attributable to the large 
uncertainties in the temperature measurements made us- 
ing Einstein, Ginga, and EXOSAT dat'd, some of it may 
also be due to unresolved substructure (e.g., cooling flows) 
and point sources. Such issues become particularly impor- 
tant when investigating the central regions of clusters as 
opposed to its integrated properties. Indeed, new high- 
resolution data obtained by Chandra support this idea 
(see, e.g., Stanford et al. 2001). We anticipate that future 
data obtained by both Chandra and XMM-Newton will 
place much tighter constraints on the Mgas{r = 0.25h^^ 
Mpc) - Tx relation and possibly even allow one to probe 
the mild break in the relationship predicted by the pre- 
heated models. 

4.3. Test 3: Mgas(r = 0.50/1"! j^p^j _ 

The solid pentagons in Figure 5 represent the gas mass 
determinations of White et al. (1997) within r = 0.50/1"^ 
Mpc using surface profile fitting of Einstein data and 
emission-weighted gas temperatures from the literature, 
for clusters with Tx > 3 keV and whose error bars are 
1 keV or smaller. Again, the predictions of the isother- 
mal self-similar model are represented by the dot-dashed 
line while the long-dashed, short-dashed, dotted, and solid 
lines represent the preheated models with Kq = 100, 200, 
300, and 427 keV cm^, respectively. 

A linear regression fit to the 20 clusters taken from 
White et al. (1997) yields a best-fit slope and intercept 
that is consistent with only the Kq = 300 and 427 keV 
cm^ preheated models (90% confidence; see Table 3). Once 
again, the isothermal model is ruled out with greater than 
99% confidence. This follows the same general trend dis- 
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covered in the previous tests. 

As expected, the influence of the gas density profile 
on the Mgas{r = 0.50/i~^ Mpc) - Tx is minimal. Only 
a very modest break is detectable at Tx « 4.5 keV 
for the Kq = 427 keV cm^ preheated model. Like the 
Mgas(r5oo) — Tx relation, this relation is mostly sensitive 
to the temperature shift (at least for clusters with Tx > 3 
keV). 

5. COMPARISON WITH PREVIOUS THEORETICAL STUDIES 

Only two other theoretical studies have examined the 
effects on the Mgas — Tx relation of entropy injection into 
the ICM: Loewenstein (2000) and Bialek et al. (2001). 
Both studies investigated the Mgas{r 500) — Tx relation and 
demonstrated that entropy injection does, indeed, steepen 
the relation, in agreement with the present work (how- 
ever, neither implemented the Mgas — Tx at fixed radii 
test). These studies suggest that models that produce an 
entropy floor with a level that is consistent with measure- 
ments of groups {Kq ~ 100 keV cm^; Ponman et al. 1999; 
Lloyd-Davies et al. 2000) are capable of matching the 
observations of even hot clusters (up to 10 keV). This is 
in apparent conflict with the results presented in §3 that 
suggest that a high entropy floor of Kq > 300 keV cm^ 
is required to match the observations of hot clusters. A 
low value of the entropy floor is also in apparent conflict 
with a number of other studies that have focused mainly 
on the Lx — Tx relation of hot clusters. For example, da 
Silva et al. (2001), Tozzi & Norman (2001) and BBLP02 
have all concluded that such low levels of entropy injection 
do not bring consistency between observations and theo- 
retical models of hot clusters. As such, a closer analysis 
of Loewenstein (2000) and Bialek et al. (2001) studies is 
warranted. 

5.1. The Loewenstein (2000) Models 

To model the observed deviations of the cluster X-ray 
scaling relations, Loewenstein (2000) has constructed a 
suite of hydrostatic polytropic models (which are normal- 
ized to observations of high-temperature clusters and nu- 
merical simulations), and then modified them by adding 
various amounts of heat per particle at the cluster center. 
Strictly speaking, the Loewenstein (2000) models cannot 
be characterized as preheated models, since the injection of 
entropy into the ICM occurs after the cluster has formed. 
Thus, a straightforward comparison between the Loewen- 
stein (2000) and BBLP02 models is not trivial. However, 
success in matching the Mgas{r^QQ)—Tx relation (the data 
of MME99) is claimed by Loewenstein (2000) for a model 
that "produces an entropy - temperature relation with the 
observed entropy floor at « 100 keV cm^." Regardless of 
how the entropy floor actually arose, this contradicts the 
results presented in §3, which suggest that an entropy floor 
of > 300 keV cm^ is required to match the observations. 
Can the analysis of Loewenstein (2000) and that of the 
present work be reconciled? 

A closer investigation of Figure 4 of Loewenstein (2000) 
reveals that flrst of all, his heated models were not com- 
pared to the actual data but rather to points that rep- 
resent MME99's best-fit power-law match to their data. 
Second, this power-law relationship was assumed to hold 
true and hence, extrapolated to span a wider range in tem- 



peratures than considered by MME99. Of the 45 clus- 
ters studied by MME999, only one had a temperature 
below 3 keV (it was 2.41 keV), yet Loewenstein (2000) 
compared his heated models to the best-fit relation of 
MME99 over the range 1 keV < Tx < 10 keV. As previ- 
ously mentioned, however, entropy injection preferentially 
affects low-temperature systems and, therefore, extrapo- 
lating scaling relations derived from high-temperature sys- 
tems down to the low-temperature regime is not safe. 

In Figure 3, we compare the best- fit heated model of 
Loewenstein (2000) (his e = 0.35 model, as the thick dot- 
long-dashed line) with the predictions of BBLP02 mod- 
els and the data of MME99. The plot clearly demon- 
strates that his best-fit model does not match the data of 
MME99 nearly as well as the BBLP02 preheated models 
with Kq > 300 keV cm^, especially at the high tempera- 
ture end. The difference in temperature ranges examined 
by Loewenstein (2000) and the present study (whose range 
of temperatures were purposely chosen to match the obser- 
vational data) has likely led to an underestimation of the 
entropy floor in these clusters by Loewenstein (2000). We 
once again re-iterate that it is extremely important that 
comparisons between theoretical models and observations 
are done over the same range in temperatures. 

5.2. The Bialek et al. (2001) simulations 

In similarity to the present work, Bialek et al. (2001) 
investigated the impact of preheating on the Mgas — Tx 
relation for a number of different levels of entropy injec- 
tion, spanning the range keV cm^ < 335 keV 
cm^. Fitting their Mgasir^m) — Tx simulation data over 
the range 2 keV < Tx ^ 9 keV, which is similar (but 
not identical) to the MME99 sample, they claim success 
in matching the observations of MME99 for models with 
entropy injection at the level of 55 keV cm^ % Kq < 140 
keV cm^, at least on the basis of slope. Their models 
with higher levels of entropy injection, apparently, pre- 
dict relations much too steep to be consistent with the 
data of MME99. These predictions are inconsistent with 
the results of the BBLP02 analytic models with similar 
levels of entropy injection (e.g., for Kq « 300 keV cm'^, 
BBLP02 predict Mgas oc T}^^'^ while Bialek et al. find 
Mgas OC T|'37) . However, we believe the difference in the 
predictions (and conclusions) of Bialek et al. (2001) and 
the present work can be reconciled. 

As noted by Neumann & Arnaud (2001), Bialek et al. 
(2001) have simulated very few hot clusters and, although 
they fit their Mgas{rnQo)—Tx simulation data over a range 
similar to MME99, the results are too heavily weighted by 
the cool clusters {Tx ^ 3 keV) to be properly compared 
with the data of MME99. As an example, we consider 
their "S6" sample of 12 clusters that have Kq — 335 kcV 
cm^. According to the present study, this model should 
give a reasonably good fit to the MME99 observational 
data, much better than that of a model with Kq « 100 
keV cm^. Although the normalization of the S6 model is 
in excellent agreement with the MME99 data (as is ap- 
parent in Table 3 of Bialek et al. and the general trends 
in their Figure 1) they rule this model out based on the 
fact that the predicted slope is 2.67, much steeper than 
the 1.98 found by MME99. However, a closer analysis re- 
veals that the fraction of cool clusters in the simulation 
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data set is much higher than fraction of cool clusters in 
the MME99 sample. For example, in the MME99 sample 
of 45 clusters, only one cluster has a temperature below 3 
keV. In the Bialek et al. (2001) S6 set, however, 5 of the 12 
clusters have temperatures below 3 keV. In addition, the 
mean temperature of clusters in the MME99 sample is ~ 
5.5 keV, while it is only about 3.8 keV in the Bialek et al. 
(2001) S6 data set. As previously mentioned, preheating 
preferentially affects low temperature systems and, there- 
fore comparisons between theory and observations should 
be done over the same range in temperatures. To illus- 
trate the problems of comparing theoretical models and 
observations that span different temperature ranges, we 
tried to reproduce the fit of Bialek et al. (2001) to their 
S6 data set. We used data presented in their Table 2 for 
clusters with Tx > 2 keV (we use their preferred "pro- 
cessed" temperatures) and fit it with a linear model and 
found Mgas oc T^-42±o.i7 rpj^jg jg slightly different from 
the value listed in their Table 3, presumably because Ta- 
ble 2 is based on data within r2oo while Table 3 is based 
on data within rsoo (they note that a change of up to 6% 
in the predicted slope can occur when switching between 
the two) . To match the conditions of the present work, we 
then discarded all simulated cluster data below 3 keV (the 
mean temperature for the remaining 7 clusters was then 
5.1 keV, similar to the MME99 data) and found a best fit 
of Mgas oc T^-99±0-30. This is excellent agreement with the 
results of MME99 and only marginally inconsistent with 
the BBLP02 models of similar entropy injection. 

What about their favored models? We have tried the 
same type of test on their S3 data set {Kq ~ 100 keV cm^). 
Fitting all simulated clusters with Tx > 2 keV (mean tem- 
perature of 3.8 keV) we find Mgas oc T^-86±o.i2^ .^j^j^^j^ jg 
in good agreement with the results of MME99. When we 
remove all clusters below 3 keV (mean temperature of 4.9 
keV), however, the best fit is Mgas oc T^-77±o.38_ 
case, the best-fit relation is not very constraining. It is 
even indistinguishable from the self-similar result. It is 
apparent from their Figure 1, however, that the predic;ted 
normalization for this model (and all other low entropy 
models) does not match the observations of MME99. This 
is noted by the authors themselves. They claim the dif- 
ference in the zero point may be resolved by reducing the 
baryon fraction by ~ 20%. As we noted earlier, however, 
a similar normalization offset is also seen in the total clus- 
ter mass - temperature (M — Tx) relation and this cannot 
be resolved by reducing the baryon fraction. This sug- 
gests that the problem lies with the temperature, rather 
than the gas mass. Alternatively, Bialek et al. also sug- 
gest that rescaling their simulations for Ho = 70 km s~^ 
Mpc~^ (instead of 80 km s~^ Mpc~^) would bring con- 
sistency between the normalization of this model and the 
observations. This would be true only if the baryon frac- 
tion was held fixed at 0.1 and not rescaled for the new 
cosmology. Given that they assume flm =0.3, this would 
imply fli, = 0.015ft-~^ which is roughly 30% lower than 
observed in r[uasar absorption spectra (Buries & Tytler 
1998). Thus, while the normalization offset between their 
theoretical model and the observations of MME99 is di- 
rectly reduced by decreasing the value of h, it is indirectly 
increased by roughly the same proportion through the in- 
creased value of Clb/Qrn- 



In summary, as with the Loewenstein (2000) models, 
we find that the difference in the results and conclusions of 
Bialek et al. (2001) and the present work can be explained 
on the basis that different temperature ranges were exam- 
ined. In particular, we have shown that the fraction of 
cool clusters in Bialek et al.'s simulated data set is much 
larger than that found in the MME99 sample and this 
has likely led to an underestimation of the entropy floor 
in these clusters. In order to safely and accurately com- 
pare the preheated models of BBLP02 with observations 
we have paid special attention to only those hot clusters 
with Tx > 3 keV. As such, we believe our comparison is 
more appropriate. 

6. DISCUSSION & CONCLUSIONS 

Motivated by a number of observational studies that 
have suggested that the Mgas — Tx relation of clusters 
of galaxies is inconsistent with the self-similar result of 
numerical simulations and by the launch of the Chandra 
and XMM-Newton satellites, which will greatly improve 
the quality of the observed Mgas — Tx relation, we have 
implemented the analytic model of BBLP02 to study the 
impact of preheating on Mgas — Tx relation. The pre- 
dictions of the model have previously been shown to be 
in very good agreement with observations (e.g., Lx — Tx 
relation and Lx — cr relation). 

In agreement with the previous theoretical studies of 
Loewenstein (2000) and Bialek et al. (2001), our analysis 
indicates that injecting the intracluster medium with en- 
tropy leads to a steeper relationship than predicted by the 
self-similar result of numerical simulations of clusters that 
evolve through the effects of gravity alone. Loewenstein 

(2000) and Bialek et al. (2001) have found that mod- 
els that produce an entropy floor of Kq ~ 100 keV cm^, 
which is consistent with measurements of galaxy groups, 
are capable of reproducing the Mgas — Tx relation of hot 
clusters. This is inconsistent with our analysis, which in- 
dicates that a "high" level of entropy injection {Kq > 300 
keV cm^) is required to match the observational data of 
hot clusters of White et al. (1997), Peres et al. (1998), 
and MME99. It is also inconsistent with BBLP02's best- 
fit value of Kq « 330 keV cm^ found via an investigation 
of the Lx — Tx relation of both groups and hot clusters. 
They note that the strongest constraints for a high entropy 
floor comes from hot chisters. Moreover, a high value of 
Kq, one that is inconsistent with the predictions of the 
best-fit models of Loewenstein (2000) and Bialek et al. 

(2001) , has also been reported by Tozzi & Norman (2001). 
Finally, da Silva et al. (2001) used numerical simulations 
with a "low" value of Kq 80 keV cm^ (which is similar 
to the predictions of the best-fit models of Loewenstein 
2000 and Bialek et al. 2001) and found that they could 
not reproduce the observed X-ray scaling relations. Our 
result, on the other hand, is consistent with the results 
of BBLP02, Tozzi & Norman (2001), and da Silva et al. 
(2001). As discussed in §5, we believe the difference be- 
tween the studies of Loewenstein (2000) and Bialek et al. 
(2001) and present work can be explained by considering 
the difference in temperature ranges studied. In particu- 
lar, we have focused only on hot clusters in an attempt 
to match the majority of the observational data as closely 
as possible. The results and conclusions of the other two 
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studies, however, are strongly influenced by their low tem- 
perature model data. 

We have proposed that the Mgas — Tx relation can be 
used as a probe of the gas density profiles of clusters if 
it is evaluated at different fixed radii. This is a new test. 
The preheated models of BBLP02 predict a mild break 
in the scaling relations when small fixed radii (such as 
r = 0.250/i~^ Mpc) are used. The scatter in the cur- 
rent observational data is consistent with the predictions 
of the BBLP02 models with Ko > 300 kcV cm^: however, 
the exact shape of the gas density profiles is not tightly 
constrained. We anticipate that large samples of clusters 
observed by Chandra and XMM-Newton will place much 
stronger constraints on the gas density profiles of clusters 
and allow for further testing of the preheating scenario. 

Finally, the high level of energy injection inferred from 
our analysis has important implications for the possible 
sources of this excess entropy. Valageas & Silk (1999), 
Balogh et al. (1999), and Wu et al. (2000) have all shown 
that galactic winds driven by supernovae can only heat the 
intracluster/intergalactic medium at the level of < 0.3—0.4 
keV per particle. This is lower than the 1 — 2 keV per par- 
ticle result foimd here. Thus, if the BBLP02 preheated 
models provide an accurate description of the ICM, su- 
pernovae winds alone cannot be responsible for the excess 
entropy. It has also been speculated that quasar outflows 
may be responsible (e.g., Valageas & Silk 1999; Nath & 
Roychowdhury 2002). This remains an open possibility. 



The role of radiative cooling also remains an open is- 
sue. Recently, it has been suggested that both radiative 
cooling and preheating together could be actively involved 
in shaping the X-ray scaling relations (e.g., Voit & Bryan 
2001; Voit et al. 2002). Radiative cooling (and subsequent 
star formation) would serve to remove the lowest entropy 
gas, which in turn would help to compress the highest en- 
tropy gas, thus increasing the emission- weighted gas tem- 
perature and steepening the Mgas — Tx relation (cf. the 
discussion of entropy in the Cool+SF simulation of Lewis 
et al. 2000). In this way, the combination of cooling and 
preheating may reduce the best-fit entropy level, perhaps 
even to a level that can be provided by supernovae winds 
(Voit et al. 2002). Further study is required to determine 
the relative roles that both preheating and cooling have 
on cluster evolution. 
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Table 1 

Results of line.vr recressiox fits to .A/,,„,s(;-,o(i) - Tx data 



Model 


Entropy Floor (keV cm^) 


m« 


b" 


Isothermal model 





1.50 


12.72 


Preheated models 


100 


1.65 


12.54 




200 


1.66 


12.48 




300 


1.67 


12.45 




427 


1.68 


12.42 


MME99 data 


? 


1.93 ±0.16 


12.21 ±0.13 


MME99 data minus 2 clusters 


? 


1.82 ±0.14 


12.31 ±0.12 


Note. — We have fit models of 


the form log Mgas = m log Tx 


± b over the 


range 3 keV <Tx < 



10 keV. 



"^Uncertainties correspond to the 90% confidence level. 



Table 2 

Results of linear regression fits to Mgas{r = 0.25ft~^ Mpc) - Tx data 



Model 


Entropy Floor (keV cm^) 


m" 




Isothermal model 





0.84 


12.61 


Preheated models 


100 


0.91 


12.43 




200 


0.95 


12.34 




300 


1.06 


12.21 




427 


1.19 


12.03 


Peres et al. and White et al. data 


? 


1.11 ±0.22 


12.06 ±0.17 



Note. — We have fit models of the form log Mgas = m log Tx ± h over the range 3 keV < Tx ^ 
10 keV. 

■^Uncertainties correspond to the 90% confidence level. 



Table 3 

Results of linear regression fits to Mgasir = 0.507i~^ Mpc) - Tx data 



Model 


Entropy Floor (keV cm^) 


m» 


b« 


Isothermal model 





0.97 


12.87 


Preheated models 


100 


0.98 


12.78 




200 


0.98 


12.75 




300 


0.99 


12.72 




427 


1.01 


12.68 


White et al. data 


? 


1.12 ±0.25 


12.53 ±0.20 



Note. — We have fit models of the form log Mgas = 'm log Tx + b over the range 3 keV < Tx < 
10 keV. 

^Uncertainties correspond to the 90% confidence level. 



